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[57] ABSTRACT 

In a new formulation for digital phase-locked loops, loop- 
filter constants are determined from loop roots that can each 
be selectively placed in the s-plane on the basis of a new set 
of parameters, each with simple and direct physical meaning 
in terms of loop noise bandwidth, root-specific decay rate, 
and root-specific damping. Loops of first to fourth order are 
treated in the continuous-update approximation (B L T — >0) 
and in a discrete-update formulation with arbitrary B L T. 
Deficiencies of the continuous-update approximation in 
large-B/T applications are avoided in the new discrete- 
update formulation. 
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METHOD OF IMPLEMENTING DIGITAL 
PHASE-LOCKED LOOPS 


ORIGIN OF INVENTION 

The invention described herein was made in the perfor- 
mance of work under a NASA contract, and is subject to the 
provisions of Public Law 96-517 (35 USC 202) in which the 
Contractor has elected not to retain title. 


TECHNICAL FIELD 


10 


The present invention comprises a method of implemen- 
tation of a digital phase-locked loop with loop constants 
derived from roots placed in the s-plane on a root-by-root 
basis in terms of root- specific damping and root specific 
decay rates. The method provides improved flexibility in 
tailoring high-order loop performance and the design is 
more straightforward and understandable for high-order 
loops. The inventive method provides a fully digital formu- 
lation, free of analog complications. Further, loop bandwidth 
and damping do not change for high-gain loops. 


BACKGROUND ART 

It is desired to simplify and improve the design, analysis 25 
and synthesis of high-order digital phase-locked loops; to 
find a technique for selecting loop constants in high-gain 
loops in a manner that leads directly to desired damping and 
loop noise bandwidth. 

30 

Previous analyses of digital phase-locked loops (DPLLs) 
are based on the traditions of analog loops and introduce 
unnecessary analog considerations such as loop-filter time 
constants and uncontrolled gain variations. This reliance on 
analog tradition makes digital-loop analysis unnecessarily 35 
cumbersome and circuitous and impedes the progress of 
analysts with little analog training. Theory for digital loops 
can be rigorously developed from first principles without 
reference to analog concepts. With an appropriately formu- 
lated “digital analysis,” one discovers that DPLL theory and 4Q 
design become more straightforward (particularly for third 
and fourth order loops) and that loop performance is more 
accurately controlled for “high gain” loops. 

Previous analyses have begun with the closed-loop equa- 
tion in the “continuous update” (CU) limit in which B^T+0, 45 
where B L is the loop noise bandwidth and T is the loop 
update interval. For sufficiently small B/T (e.g., B^T^iO.02), 
the CU approximation can provide an adequate starting 
point for loop analysis and design. When B^T is increased 
in this model to larger, “high-gain” values, however, loop 50 
roots can move away from their initial “small B L T’ paths in 
unplanned directions and the loop can diverge from 
expected behavior. For a loop with discrete update (DU) 
intervals, a solution to the loop equation can be developed 
in which root locations follow predetermined paths as a 55 
function of B L T. This feature, which is an automatic benefit 
of a new controlled-root parameterization of the present 
invention, can provide, for example, supercritically-damped 
response for all allowed B L T values. 

The following U.S. patent numbers were found in a search 60 
of relevant prior art but none is deemed to affect the 
patentability of the inventive method hereof. 

U.S. Pat. No. 3,772,600 Natali 

U.S. Pat. No. 4,426,712 Gorski-Popiel 65 

U.S. Pat. No. 4,769,816 Hochstadt et al 
U.S. Pat. No. 4,771,250 Statman et al 
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U.S. Pat. No. 4,794,341 Bartonetal 
U.S. Pat. No. 4,820,993 Cohen et al 
U.S. Pat. No. 4,821,294 Thomas, Jr. 
U.S. Pat. No. 4,890,305 Devries 
U.S. Pat. No. 4,975,930 Shaw 
U.S. Pat. No. 5,016,005 Shaw et al 
U.S. Pat. No. 5,036,296 Yoshida 
U.S. Pat. No. 5,073,907 Thomas, Jr. 
U.S. Pat. No. 5,109,394 Hjerpe et al 
U.S. Pat. No. 5,122,761 Wischermann 


STATEMENT OF THE INVENTION 

The instant invention comprises an approach that reverses 
the conventional procedure of design for obtaining loop 
filter constants. The instant innovation is a method in which 
loop roots are first placed in the s-plane on the basis of new 
root-specific damping and root- specific-decay -rate param- 
eters. Loop constants are then calculated on the basis of 
these roots. 

In the present invention, loop-filter constants are specified 
in terms of new loop parameters. By design, each of these 
parameters has a simple and direct physical meaning in 
terms of a useful loop property: loop noise bandwidth, 
transient decay time or damping. For example, a simple 
choice of parameter values will automatically give a loop a 
particular selected loop bandwidth and supercritically- 
damped behavior (i.e., all roots real, negative, and equal). 
Thus, the need to solve for root location as a function of 
standard loop parameters (e.g., B^, r and k for a third-order 
loop) is eliminated and analysis is simplified. The new 
parameterization is made feasible in a practical sense by the 
fact that digital loops can usually be designed so that they do 
not suffer significantly from the effects of gain variations. 
That is, variations in signal amplitude, due to either gain 
instability or signal-power changes, can usually be 
accounted for by using a normalized phase extractor. So 
comprised, a “Hilly digital” DPLL does not require the 
analysis or precautions necessitated in other DPLL designs 
by potential gain variations. A particularly appealing benefit 
of ffie invention is the elimination of the equivocal practice 
of using the loop parameter r as both “damping factor” and 
an overall gain factor. 

The analysis is extended to fourth-order loops because of 
the potential advantages of such loops. In some spacecraft 
applications, loop bandwidth can be set to a smaller value 
for a fourth-order loop than for a third order loop. Conse- 
quently, lower signal strengths can be reliably tracked. 
Fourth-order DPLLs, unlike their analog counterparts, are 
easy to design and implement, given the new fully-digital 
formulation. Accurate placement of loop roots results from 
a simple selection of parameters values rather than compli- 
cated analog circuit design. 

To establish a foundation for analysis, a high-level 
description of a DPLL is presented by way of background. 
For loops of first to fourth order, the new controlled-root 
parameterization is used to derive a CU-limit solution and to 
develop a general approach from which numerical, con- 
trolled-root solutions to the DU loop can be derived. Solu- 
tions are given for phase and phase-rate feedback, with 
computation delay for closing the loop set to either zero or 
one update interval. To tie in with traditional analysis, the 
new loop parameters are related to old loop parameters in the 
CU limit. Loop transfer functions are presented for each 
loop order. The technique for direct transient-free acquisi- 
tion is extended to fourth-order loops. Finally, results for two 
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measures of loop performance, mean time to first cycle slip 
and steady- state phase error, are described. 

OBJECTS OF THE INVENTION 

It is therefore a principal object of the present invention 5 
to provide a method for specifying loop- filter constants of 
digital phase-locked loops in terms of loop parameters 
having direct physical meaning as a useful loop property. 

It is another object of the invention to provide a method 
of simplifying the analysis of digital phase-locked loops by 10 
obviating the need to solve for root location as a function of 
standard loop parameters. 

It is still an additional object of the present invention to 
provide a method for parameterizing digital phase-locked 15 
loops, wherein loop constants are computed on the basis of 
loop noise bandwidth and controlled-root independent 
parameters related to decay times and damping. 

BRIEF DESCRIPTION OF THE DRAWINGS 

20 

The aforementioned objects and advantages of the present 
invention, as well as additional objects and advantages 
thereof, will be more fully understood hereinafter as a result 
of a detailed description of a preferred embodiment when 
taken in conjunction with the following drawings in which: 2 5 

FIG. 1 is a block diagram of a digital phase-locked loop 
with phase and phase-rate updates; 

FIG. 2 is a graphical illustration of the relation between 
the conventional second-order loop filter parameter r and the 
controlled-root parameter rj 1 1 in the continuous update limit; 30 

FIG. 3 is a graphical illustration of relations between 
third-order loop filter parameters r, k and controlled-root 
parameters ^ in the continuous update limit with the 
shaded area representing values of r, k where all roots are 
real; 35 

FIG. 4 is a graphical illustration of the relationship 
between fourth-order loop filter parameters r, K, a and 
controlled-root parameters r| 2 , t\ 2 in the continuous 

update limit; 40 

FIG. 5 is a graphical illustration of transfer functions for 
loops of order 1 to 4 for various values of T] 2 and 

FIG. 6, comprising FIGS. 6a and 6b , is a graphical 
illustration of normalized loop bandwidth B L T versus ref- 
erence decay rate parameter B L T for a third-order loop 45 
without, and with computation delay, respectively; 

FIG. 7, comprising FIGS, la and 7b, is a graphical 
illustration of transfer functions for first and second order 
loops with phase and phase-rate feedback and no computa- 
tion delay; 50 

FIG. 8, comprising FIGS. 8 a and 8 b, is a graphical 
illustration of transfer functions for third and fourth order 
loops with phase and phase-rate feedback and no computa- 
tion delay; 

55 

FIG. 9, comprising FIGS. 9 a and 9b, is a graphical 
illustration of transfer functions for first and second order 
loops with phase and phase-rate feedback and computation 
delay of one update interval; 

FIG. 10, comprising FIGS. 10a and 10 b, is a graphical 6 o 
illustration of transfer functions for third and fourth order 
loops with phase and phase-rate feedback and computation 
delay of one update interval; 

FIG. 11, comprising FIGS. 11a and lib, is a graphical 
illustration of root-locus plots for a second order loop using 65 
discrete-update formulation (thick line) and continuous- 
update approximation (thin line) as a function of B^T; 


4 

FIG. 12, comprising FIGS. 12a and 12b, is a graphical 
illustration of mean time to first cycle slip for discrete- 
update loops with supercritical damping, phase and phase- 
rate feedback, no computation delay and sine phase extrac- 
tor; 

FIG. 13, comprising FIGS. 13a and 13b, is a graphical 
illustration of quantities describing steady-state phase error 
in discrete-update loops with phase and phase-rate feedback 
and no computation delay; and 

FIG. 14, comprising FIGS. 14a and 14b, is a graphical 
illustration of quantities describing steady-state phase error 
in discrete-update loops with phase and phase-rate feedback 
and computation delay of one update interval. 

1.0 DETAILED DESCRIPTION OF THE 
INVENTION 

2.0 HIGH-LEVEL DESCRIPTION OF A DPLL 

The block diagram in FIG. 1 shows the basic elements of 
a DPLL. The example DPLL shown in FIG. 1 is based on 
“immediate update” of the loop filter and feedback of phase 
as well as phase-rate. Alternate DPLL designs might feed- 
back only phase rate and/or have a substantial computation 
delay (“transport lag”). The analysis below also treats the 
case when the computation delay is equal to one update 
interval. 

An incoming signal is sampled at a high rate ( f s ) and then 
counter-rotated sample by sample with model phase values 
generated by a number-controlled oscillator (NCO) as 
directed by loop feedback. The resulting counter-rotated 
signal, which should have very low frequency, is then 
accumulated over an update interval length T. A phase 
extractor then processes the resulting complex sum to pro- 
duce a value for residual phase for the given interval (nth). 
For an ideal phase extractor, the n lh residual phase is 
equivalent to the difference of the n th input signal phase (j)„ 
and the n th model phase (j> m n : 

( 2 . 1 ) 

with each reference to the center of the sum interval. Even 
though actual residual phase can deviate from this linear 
model due to nonlinearity in the phase extractor and/or cycle 
ambiguities, the theory presented here will be based on the 
approximation of Equation (2.1). This n th residual phase is 
immediately passed to the loop filter to assist in the predic- 
tion of the phase rate for the (n+\) th interval. The loop filter 
generates an estimate of phase rate, (j) mn+1 T, in the form of 
phase change per update interval. An estimate of the next 
model phase, the (n+iy h , is projected ahead to interval 
center by adding this estimate phase change to the previous 
n th model phase: 

( 2 . 2 ) 

The (n+iy* model phase, along with estimated phase rate, is 
used to initialize the phase and rate registers of the NCO in 
a manner that causes the NCO to produce over the (n+l) r;z 
interval, phase values characterized by said rate and center- 
interval phase. In this manner, the loop is closed and a new 
value for residual phase is produced to repeat the process. 

3.0 LOOP THEORY 

3.1 LOOP FILTER 

As suggested above, residual phase from a given interval 
can either immediately enter the loop filter and influence 
feedback in the very next interval or its effect can be delayed 
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to allow more time for computation. Two cases are consid- 
ered here: a) “immediate” updates with acceptably small 
computation times and b) updates delayed by one update 
interval. Immediate updates are discussed first. 

The use of residual phase from a given interval to update 5 
NCO phase in the very next interval can lead to a “dead 
time” during which necessary update computations are 
performed, and during which no counter-rotated phasors are 
accumulated. If the length of this “dead time” can be made 
such a small fraction of an update interval that the SNR loss 10 
is acceptably small, the immediate update of the loop filter 
becomes feasible. For immediate updates, an N th -order 
digital loop filter, using residual phase values 5<|> f up to and 
including the n th sum interval, estimates phase rate for the 
(n+iy* interval according to 15 
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-continued 

n i nii 

+ K 3 I I X T&}>/ + K 4 I T X T X Tbtyk + . 
i=l j = 1 r=l 7=1 h=l 


where estimated phase rate has been rewritten on the basis 
of Equation (2.2) under the assumption that the NCO is 
updated in both phase and phase rate. In the limit T— >0, the 
first term becomes a derivative and the sums become inte- 
grals: 


<f> m = Ki5<J) + K 2 


f 1 <//5<|) + K 3 f / [ 1 i 
J to J t'J t 0 


•fJT 

J t J t J to 


n,n+\T ~ jRTiS^n + K’l X b$i + K s X 2 8<(>7 + 
i=l i=l 7=1 


£4 X X X % + ... 20 
i=l 7=1 k=l 


where <j> m , n+1 T is phase change per update interval T and 
the loop constants K t are to be specified below. (An N'* 
-order loop has N terms.) 

In the case of a relatively long computation time, loss of 
data can often be eliminated by performing filter computa- 
tions in parallel with signal accumulation in each update 
interval, and applying the result one interval later than the 
immediate-update approach. For a computation delay 
(“transport lag”) of one update interval, the loop filter 
estimates phase rate by computing during the n th interval the 
expression 


n - 1 
1=1 


n -1 i 

h + Ks X X 5<j >;H 

f=l 7=1 


n~ 1 i i 

K a £ I X 6 <bjt h 

i=l 7=1 fc=l 


Note that the only change between Equations (3.1) and 
(3.2) is that the n th residual phase is not used in computing 
the (n+l) /; * phase-rate estimate. 

In the following discussion, the loop equation will be 
solved on the basis of immediate updates Equation (3.1), 
while only results are presented for Equation (3.2). 


3.2 CONTINUOUS-UPDATE LIMIT 


3.2.1 LOOP EQUATION 

In many applications, the update interval T is much 
shorter than all other filter time scales, and considerable 55 
insight may be gained by writing Equation (3.1) in the 
“continuous-update limit”, T^O. (One can easily show that 
Equation (3.2) leads to the same result in this limit.) To 
facilitate this, define new loop constants K f as 

60 

K= ki T for i=l, N (3.3) 

so that Equation (3.1) becomes 




^=Ki5(j) rt + K2 £ 7 tyi + 
1 i = 1 


Thus, in this limit, the basic equation governing “NCO 
rate” is the same equation that governs the VCO rate of an 
analog loop, given perfect integrators. 

Solutions of the closed-loop equation can be used on the 
theory of differential equations after substituting Equation 
(2.1) and differentiating N-1 times: 


dN * -v dN ~ l 

dt N dt N ~ l 


„+. . .+KN$ m 


<l> + *2 <|> + . . . + Kjv<j> 

dr 1 


Solution of this differential equation will give the behavior 
of model phase in response to input phase $. 

3.2.2 TRANSFER FUNCTION AND LOOP 
NOISE BANDWIDTH 

To find the frequency response of the loop, take the 
Fourier transform of both sides of Equation (3.6), and utilize 
the relation 


f d^(t) 1 

F j > = ( f 27cv )”/ r { 4>(0 } 


where F {} represents a Fourier transform. This produces 

[(i27rv) N +K 1 (i2jcv) yv “ 1 +K 2 (i2jcv) w_2 + . . . +K JV](( , m (v)=lK 1 (i2jcv) Ar - 1 + 
K 2 (i27EK) W ~ 2 + . . . +K^](j)(v) (3.S) 

where 0(v) and <j) m (v) are the Fourier transforms of <j>(t) and 
(t). The closed-loop transfer function H(v) is defined by 
<j) m (v)=H(v)<j)(v), so we have 

Ki(i27cv)^ _1 -(- K2(i27rv)^ 2 + . . . -hcjv (3-9) 

H(y) = 

(Hkv) n + Ki(i2jcvy^ -1 + k 2 (i‘ 2 jiv ) 2V ' _2 + . . . + k n 

Plots of the transfer function for various loop orders and 
constants are shown in FIGS. 7 to 10. The single-sided loop 
noise bandwidth for the closed loop is defined by 


25z,sJ dv\H(y)} 2 

This integral is evaluated in Appendix A to find B L as a 
function of xq, k 2 , . . . continuous-update limit and the 
results are summarized in Table 3-1 for loop fourth order. 

3.2.3 PARAMETERIZATION OF LOOP 
CONSTANTS 

Loop bandwidth provides us with one parameter for 
parameterizing the loop constants. Conventional loop theory 
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specifies additional parameters at this point (e.g., r, k, . . . ) 
to complete the parameterization. However, since r, K, etc. 
can not be given direct physical significance for loops higher 
than second order, they are less useful than more carefully 
chosen parameters. 5 

To begin the process of defining new parameters, define 
the higher k’s in terms of iq: 


KpOt/Kj f for i=2,N (3.11) 

No loss of generality is suffered at this step since three new 10 
parameters replace the three k’s. When this equation is 
substituted for the k’s in the expression for loop bandwidth, 
one obtains as a function of the a’s, as shown in Table 
3-L This equation can be solved for K a in terms of B^ and 
the a’s, as shown by the first line of equations for each loop 
order in Table 3-2. Based on this expression for K* and 
Equation (3.11), the loop constants can be expressed as a 
function of B^ and the a’s. 

The variables a 2 , . . . , a^, will now be parameterized in ^ 
terms of physically significant quantities. To begin the 
search for such parameters, return to Equation (3.6). When 
the roots are unequal (i.e., non- degenerate), the general 
solution to the homogeneous equation (<()(t)=0) is of the form 

N (3.12) 25 

<W0 = L 

i=l 

where each s f is a root of the characteristic equation 

S^+Kj 1 -hc 2 s w-2 + . . . +K,v=0 (3.13) 

30 

and where the a f are amplitudes to be determined by the 
initial conditions. We can use characteristics of these roots 
to parameterize k 1s k 2 , . . . , k^, so as to create a more 
physically interpretable representation. 

Since the loop constants k x are real, we can parameterize 35 
the roots without loss of generality as 


{Si,s 2 ,s 3 j s 4 ,s s ,s 6 ; . . . M-PjUtru); - 

p 2 (i±n 2 ; -p 3 (i±n 3 ); . . . } (3.i4 

9 n 40 

where q ? - and p f * are real. For an odd number of roots, the 
unpaired root is not given an r| 2 parameter. This parameter- 
ization takes advantage of the fact that complex roots of a 
polynomial equation always occur in conjugate pairs when 
the polynomial coefficients are real. Our goal here is to 45 
create parameters which dictate the relative behavior of the 
roots once B^has been specified. Thus we can choose one 
root factor, p 1# as the reference decay-rate parameter, which 
is to be determined below as a function of B^, and form new 
parameters which indicate magnitude relative to the first: 5Q 

P=X,p 1 fori^2 (3.15) 

Our root parameterization is now given by: 

{Si, s 2 ; s 3 , s 4 ;s 5 , s 6 ; . . . }={— Pi(liHi); 55 

-p^d+Ha); . . . } (3.16) 

The overall scale factor p 1? which we choose to be a positive 
real number for loop stability, is proportional to B L in the CU 
limit, as can be inferred from the specific example in Section 6Q 
3.2.4. 

The q f 2 , which we will refer to as the damping param- 
eters, contain useful information about each root pair. The 
sign of q t 2 for a given root pair determines the damping 
behavior for that root pair: 65 

q^O two real roots; overdamped 
q,. 2 =0 two real, equal roots: critically dampened 
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q 2 <0 complex conjugate pair: underdamped/oscillatory 
and the magnitude of q* 2 is a measure of the separation of 
the root pair: 



The \ ’s, which we will refer to as the relative magnitude 
parameters, control the relative magnitudes of different root 
pairs. Furthermore, a negative (positive) X 1 for a given root 
indicates an exponentially growing (decaying) solution to 
(3.6), since p 1 is chosen to be a positive real value. Some 
interesting values of A,’s and q’s are: 
all 1 all real parts of roots equal 
all q 2 =0, \=l all roots real and equal (supercritically 
damped) 

all q 2 =— 1, X=1 “standard” underdamped response 
where “standard” underdamped response for a given root 
corresponds to the response of a 2 tld -order loop with 
r=2(£=0.707). 

To express the k’s as a function of these new root 
parameters, equate each term in Equation (3.13) with the 
like term in the same polynomial factored into its roots: 


(S-S!)(S-S 2 ) . . . (s-s„)=0 

(3.18) 

k’s are then given by 


N 

(3.19) 

K\ =-l,Si 
i 


N 

(3.20) 

K2 = + Z SiSj 


i<j 


N 

(3.21) 


k 3 = - I SiSjSk 
i<j<k 


k n = (-1 ?*(5jS 2S 3 . . . s N ) (3.22) 

When the root expressions in Equation (3.16) and the K 
expressions in Equation (3.11) are substituted into these 
equations for the k’s, one can solve for p 1 and a’s in terms 
of B^, and q 2 . Results for the a’s are given for loops of 
order 1 to 4 in the second line of equations for each loop 
order in Table 3-2. To tie in with traditional parameters, the 
a’s are al so given in terms of r, q, and a in the third line in 
Table 3-2. A graphical mapping between the new parameters 
and traditional parameters is presented in FIGS. 2, 3, and 4 
for second-, third- and fourth-order loops, respectively. For 
reference, plots of the transfer functions for loops of order 
1 to 4 are presented in FIG. 5 for various values of 
and q ( - 2 

Since the transient response of a loop is characterized by 
the solution. Equation (3.12), to the homogeneous equation, 
knowledge of root locations provides a basis for predicting 
transient behavior and settling time. Because the q 2 and X t 
values, along with the loop bandwidth, completely specify 
the roots by location in the complex plane, loop transient 
response is directly selected at the outset when the new loop 
parameters are chosen. For loops of first to fourth orders, 
Table 3 — 3 presents loop constants for two likely imple- 
mentations: a) standard underdamping, where all roots have 
the same decay time (all ^=1) and all q f 2 =-l and b) 
supercritical damping, where all roots have the same decay 
time and are critically damped (all q 2 =0). For comparison. 
Table 3 — 3 also presents the corresponding conventional 
parameters. 
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3.2.4 DETAILED DERIVATION FOR A 
FOURTH-ORDER DPLL 
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3.3 LOOPS WITH DISCRETE UPDATE 
INTERVALS 


To illustrate the method hereof, a fourth- order loop will be 
treated in detail. As indicated by Equation (3.16), the roots 5 
for a fourth-order loop are expressed in terms of the con- 
trolled-root parameters as 

{Sj, S 2 ; ±IU); -JWl^E)} (3-23) 

10 

Using Equations (3.19) through (3.22) and Equation (3.23), 
we can write: 


Ki = 

-(Si + 5 2 + 53 + 54) 

(3.24) 

= 

(2 + 2 tc 2 )Pi 

15 

k 2 = 

5]S2 + 5 1 5 3 + S 1 S 4 + S2S3 + 5254 + S 3 S 4 

(4n 2 + (l-Ti, 2 ) + V(l-n 2 2 ))Pi 2 

(3.25) 

*3 = 

~(5l5 2 53 + 5l5 2 54 + Sl.S 3 .S 4 + 52 S 3 S 4 ) 

(2X2(1 -T|l 2 ) + 2X 2 2 (1 — Hn 2 ))Pl 3 

(3.26) 20 

K 4 = 

5152^3^4 

X 2 2 (l - T|i 2 )(1 - T) 2 2 )Pl 4 

(3.27) 


Equation (3.24) can be used to express in terms of iq: 
o K 1 (3.28) 

This expression and Equation (3.11) can be substituted in 
Equations (3.25) to (3.27) and the three resulting equations 
can be solved for the cc’s to give 


4>l2 + (1 - T) 1 1 2 ) + 7l2 2 (l ~ T]2 2 ) 
(2 + 2 X 2) 2 


(3.29) 

35 


3.3.1 LOOP EQUATION 

A loop with phase and phase-rate feedback and immediate 
updates will be analyzed in the text while Appendix B 
summarizes differences found in the analysis of loops with 
a computation delay of one update interval. Equation (2.1) 
can be used to rewrite Equation (3.1) in terms of the model 
phase and signal phase: 

n n i (3.34) 

A<t> m , n+ l + + #2 , £ <j) mj j + K 3 Z X § m j + . . . = 

f=l 1=1 /= 1 

n n i 

+ Kj Z (j)i + J^3 X Z 0; + ... 
1=1 i~ 1 /= 1 

where the operator A is defined by 

Ax„-X„~X n .i (3.35) 

and where Equation (2.2) has been used to replace estimated 
phase rate <j) mn+1 T with A(Jj m n+1 , under the assumption that 
the NCO is updated with both phase and phase rate. To 
convert Equation (3.34) to a difference equation, apply the 
A operator (N-l) times: 

A^1+K^ WK Z A J ^+ . . . +K^ m , w =K 1 A w -> n + 
K 2 A Ar ~ 2 <j)„+ . . . +K^ n (3.36) 

where N is the loop order. In analogy with Equation (3.6), 
solution of this difference equation will give the behavior of 
model phase (j) m in response to signal phase <j). 

3.3.2 TRANSFER FUNCTION 


2X 2 (i-ni 2 ) + 2^ 2 (i-n2 2 ) 
(2 + 2 A . 2) 3 


(3.30) To find the frequency response of the loop, take the 

z-transform of both sides of Equation (3.36) to obtain 


a31 > 

OtA 

(2 + 2X 2 ) 4 

Analytical evaluation of the loop bandwidth integral in 
Equation (3.10) gives loop bandwidth in terms of the k’s: 


Ki 2 K 2 K3 - K1K3 2 - Ki 3 K 4 + 
K 2 2 JC3 ~ XTiJCiJC* - JC 3 K 4 
4(KiK 2 K3-K3 2 -Ki 2 K4) 



Substitution of Equation (3.11) into this expression allows a 
solution for iq in terms of B L and the cc’s: 


OC 2 Ot3 -OL3 z -(Xa 


O 2 O 3 - - 04 + ct2 0C3 - 0C2C<4 - C^CX* 


. . . +K^ m (zMl-z- 1 )"- 1 K 1 

(f)(z)+ . . . +K*#(z) (3.37) 

where (j) m and <j) are the z-transforms of and (j) n , 
respectively. To reach this expression, we have used the 
relations 

Z{ & N %„ }-( 1 — z _ 1 ) N Z{%„ } (3.38) 

and 

Z{X« + il=z Z {*„} (3.39) 

where Z{} represents a z-transform. Since the closed-loop 
transfer function is defined by 


Note that this equation in combination with Equation (3.28) 
shows that p! is proportional to B L . 55 

Thus, if r^ 2 , ^ 2 , and r | 2 2 are specified, one can calculate 
the a’s using Equations (3.29) to (3.31). These values of cq 
can then be substituted in Equation (3.33) along with B L to 
obtain a value for iq. This value for iq , along with the a’s, 
lead to values for k 2 , k 3 , and k 4 through Equation (3.1 1 ). As 60 
shown in Equation (3.3) and listed in Table 3-2, the k’s are 
multiplied by an appropriate power of T to obtain the k’s to 
be used in the loop filter. In this manner, the loop constants 
for a fourth-order loop (in the limit T->0) can be calculated 55 
in terms of loop bandwidth and the con troll ed-root param- 
eters t\ 2 and 



The frequency response of the loop is obtained by substi- 
tuting 

z=e a ™ T (3.43) 
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-continued 


11 

in Equation (3.41) where v is the frequency in H z , Plots of 
the transfer function for various loop implementations are 
presented in FIGS. 7-10. 


3.3.3 LOOP-NOISE BANDWIDTH 

In analogy with Equation (3.10), the single-sided noise 
bandwidth for the closed loop is defined by 


2 B l = 


C 


T 

i ltf(e' 2nv7 )l 2 rfv 

"T 


(3.44) 10 


which can be rewritten as a contour integral in the form 


2B l t=^- § H(z)H(r l )r l dz 


(3.45) 


15 


where the closed path is along the unit circle. This integral 
can be computed on the basis of residues within the unit 
circle to obtain B L T as a function of the poles of the 20 
integrand. For simple poles, the integral is given by 


2 b l t=z {{ Z - zi )ff(z)ff(r ] )r l W (346) 

i 

where the z t are poles of the integrand within the unit circle. 2 s 
(For cases with poles of order greater than first, the residue 
evaluation must be appropriately modified.) 

As seen in Equation (3.41), the poles of H(z) are the roots 
of the polynomial D(z) in the denominator of H(z) and 
therefore satisfy the characteristic equation 30 


i K^O 


K 3 +... +z?~ 

(3.47) 


Let the N roots of this polynomial be z r - so that this equation 
can also be written as 


D(z)=(z-z x){z'Z^(z-z 3 ) . . . (z-z N )=0 (3.48) 

Since the N roots of D(z) must all lie within the unit circle 
if the loop is to be stable (see next subsection), all poles of 40 
H(z -1 ) will lie outside the unit circle and will not contribute 
to the contour integral. Furthermore, one can easily show 
that the residue for the pole at z=0 is zero since H(z~ ! )^0 
as z-^0. Thus the contour integral can be evaluated in a 
straightforward though algebraically tedious fashion on the 45 
basis of the N poles of H(z). The resulting expressions for 
B l T as a function of the roots z t are lengthy and uninfor- 
mative, particularly for the higher-order loops, and will not 
be presented here. 

50 

3.3.4 LOOP CONSTANTS FROM LOOP ROOTS 

To obtain the relationship between the roots and the loop 
constants, first collect terms according to the power of z in 
Equations (3.47) and (3.48). When the coefficients of like 55 
powers of z are equated, one obtains N equations relating 
roots and loop constants: 


1 

1 

£ 

1 

1 

- 

II 

(3.49) 


(3.50) 

. - Ajv_1 

)»- 

(3.51) 


( N ~2 )*— -*w 


ziZz^zn _( N \ ( N- 1 ( 3 - 52 ) 

z i Zi (a'-I ) (^-2 ) Kl 2 

(3.53) 

It Zi=l-Kl 
I 

where (N/k) is the binomial coefficient. These N equations 
can be used to solve for each of the N loop constants in terms 
of the N roots, z { . First solve Equation (3.53) for K x . The 
result can then be substituted in Equation (3.52) to allow a 
solution for K 2 in terms of the z f . Proceeding sequentially in 
this manner through the rest of the equation, one can obtain 
an expression for each loop constant in term of the roots. 
Thus, if the roots are known, the loop constants can be 
calculated. 

When the contour integral for B^T is evaluated as a 
function of roots for a given loop order, the result can be 
reduced to a form that contains only the functions of z £ found 
on the left-hand side of Equations (3.49) to (3.53). When this 
form is reached, B L T can then be easily expressed as a 
function of only the loop constants. The results are presented 
in Table 3-4 for loops of order one to four. 


3.3.5 SOLUTION TO HOMOGENEOUS 
EQUATION 

Solutions to the homogeneous form of Equation (3.36) 
provide information as to the transient behavior and stability 
of the loop. The general non- degenerate solution to the 
homogeneous equation is given by 


N 

§mji ~ ^ QiZt 1 


(3-54) 


where the N amplitudes are to be determined from initial 
conditions and where the N complex numbers z f are again 
the roots of the polynomial D(z) in Equation (3.47). (To see 
that these roots provide solutions to the homogeneous equa- 
tion, substitute § mtn =z n in the left-hand side of Equation 
(3.36), with the right-side set equal to zero, and reduce to the 
form of Equation (3.47).) Thus, the roots from the homo- 
geneous equation are also the poles of the transfer function. 
In order for the loop to be stable, the loop constants K f must 
be set to values that cause all the roots to fall within the unit 
circle. With a modulus less than 1, a root cannot cause the 
homogeneous solution in Equation (3.54) to diverge. 


3.3.6 LOOP PARAMETERIZATION 

Loop parameterization in the case of discrete update 
intervals parallels Section 3.2.3 for the CU limit. Loop noise 
bandwidth B^ and the same root-location parameters are 
adopted as independent loop parameters. The roots will be 
parameterized in the form 

{z lf Z 2 > z 3 , z 4 ; . . . y={e -Pl(1±T11)r ; e _X2 ^ l(1± P 2)T ; . . . } (3.55) 

where and rj £ are the N-l independent parameters speci- 
fied in Section 3.2.3, with ^=1. These parameters and 
“normalized” loop bandwidth B^T will comprise the N 
independent loop parameters needed to completely specify 
the loop. 



5,602,883 


13 14 


In Equation (3.55), the reference decay-rate parameter p is 
which will be represented in the normalized, dimensionless 
form p x T in the DU case, must be determined as a function 
of these N parameters. In the CU limit, determination of 
in terms of B L , and T| f , could be carried out explicitly, as 5 
shown in Section 3.2.3. For DU loops, however, the com- 
plexity of the equations makes a closed-form solution in the 
general case impractical. Thus, a numerical solution has 
been carried out by first selecting a value for and the 
N-l independent parameters, and r\ f , and then computing to 
numerical values for the N roots z-through use of Equation 
(3.55). The resulting z f values can be used to compute the 
normalized loop bandwidth, B/T, as shown in Section 3.3.3, 
and the loop constants, as shown in Section 3.3.4. Repeating 
the process in this fashion on the basis of the same and r^. 15 
values, one can vary the parameter p x T numerically to 
obtain B/T and the loop constants as a function of pjT. 

In general, B 7 T increases as p T T increases from zero but 
can go no higher than a loop-specific maximum value. Plots 
of B^T versus pjT are shown in FIGS. 6a and 6b for two 20 
supercritically-damped third-order loops with phase and 
phase-rate update, one with a computation delay of zero, the 
other with a computation delay equal to one update interval. 

In the zero-computation delay case, B/T can get no higher 
than 9.5, which corresponds to a p/T value of +«>. in the 25 
other case, B^T reaches a peak at PjT=-ln(%) as P X T 
increases, where B^T^O.3 is the maximum attainable value. 

For a given B L T, therefore, one can find the corresponding 
PiT, if any, and the corresponding loop constants on the 
basis of the analysis outlined above. Thus, loop constants 30 
can be determined for given X g and as a function of B^T. 
Results are presented in Tables 3-5 and 3-6 for loops of 
order one to four, given phase and phase-rate feedback, 
supercritical damping and standard underdamping. Once 
and the loop roots are known, the transfer function can be 35 
computed on the basis of Equation (3.41). Results are 
plotted in FIGS, la and lb through FIGS. 10a and 10Z? for 
the two standard loop configurations. 

Even though a general solution to the DU loop has not 
been obtained, the equations can be expanded as a power 40 
series in B L T to obtain the higher-order terms relative to the 
CU limit. Results are presented in Table 3-7 and 3-8 for the 
loop constants as a power series in B^T for the standard 
loops. 


constants in Table 3 — 3.) Where the curves separate, the CU 
approximation starts to diverge from standard underdamp- 
ing and loop transient response changes. The inset plot 
illustrates this divergence in terms of the corresponding 
damping parameter r| 2 , which starts at the intended under- 
damped value of -1 at 6^=0, increases to a critically- 
damped value of 0 at true B L T=0.8, and then approaches+1 
at true B L T~1 .2. Thus, loop response at high B^T values will 
not match original intended response. In contrast, as indi- 
cated by the thick lines, the DU exact solution maintains 
standard underdamping (i.e., rj ■ 2 =-l) for allowed values of 
B l T. 


4.0 TRANSIENT-FREE ACQUISITION WITH 
DPLLS 

If the signal phase and its time derivatives are accurately 
known at start-up, it is possible to initialize the loop sums 
and loop phase so that the loop starts tracking in-lock, with 
no transients. In steady-state tracking, residual phase 
becomes a constant so that Equation (3.1) or 

Equation (3.2) becomes 

n~n c (4.1) 

A$ n+l = K^ ss + K 2 X 50 ( + 
i=l 

n-n c i n-n c ; 

K 3 X X 5d); + /r 4 X X X % 
i= 1 j= 1 7 i=l 7=1 k=l 

for up to a fourth-order loop. The quantity n c is the com- 
putation delay, which is zero or one update interval for the 
present analysis. Note that estimated model phase rate, 
§m,n+ iT h as been replaced by differenced input phase since 
model phase tracks input phase exactly, except for a constant 


offset, in steady- state tracking. Based on this 
higher-order differences become 

expression, 

n-n c n-n c ( 

A 2 (J> n+ i = K 2 ty S s + K 3 X &j> f + K a X X 6 % 
i=l f=l 7=1 

(4.2) 

n-n c 

A% + i = K 3 fy ss + K A X Stfr 
r=l 

(4.3) 

A\ + i = 

(4.4) 


3.4 EXAMPLE OF STRAYING ROOTS AT 
LARGE B L T VALUES IN THE 
CONTINUOUS-UPDATE APPROXIMATION 

As B l T increases in the CU approximation for loop 50 
constants (e.g., Table 3-3), the path of the loop roots can 
stray from their original intended course and true loop noise 
bandwidth can exceed the B^ parameter value used to 
compute the loop constants. This deviation is illustrated in 
both the sT-plane and z-plane in FIGS. 11a and 11 b for a 55 
second-order DU loop with phase and phase-rate feedback, 
standard underdamping, and no computation delay. The 
thick straight lines, which are based on the exact DU 
solution presented in preceding subsections, show the paths 
taken by standard-underdamped roots as (true) B^T 60 
increases from 0 to 1.4. The thin curved lines, which are 
produced by the CU approximation, follow the DU lines 
until about B L T=0.1 and then curve toward the real axis. 
(Note that the CU-approximation curves are marked by both 
true B l T values and “parameter” B^T values. True B L T is 65 
the actual effective noise bandwidth for a DU loop, while 
“parameter” B L T is the value used to compute the loop 


One can easily show the differences 
phase values by 

are related to input 


(4.5) 


(4.6) 

^N + l-^n+^n~l~^n~2 

(4.7) 

A4{ t , n+l - 4<{> rt +6({) n _ 1 -4(j)„_ 2 -K|)„_3 

(4.8) 


where each phase value is referenced to interval center. We 
can relate these phase values to derivatives of incoming 
signal phase <|) Cl) by means of a Taylor expansion whose 
origin is the center of the n th interval (t=t n ): 

m (t-tn) 2 m (4-9) 

<K0 = 0/i + 6 - tn)$ } + — + 


This expansion shows that phase at the center of the nearby 
intervals is given by 



15 
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4>n+l = Wn + T) = <j)„ 4- T(J) ( n l) 


XL 

2 




4>n = <KO = <l>n 

4>«- 1 = tfe. - 31 = <t>« - ^ 2) - 


(4.10) 


XL 

6 


(4.11) 

(4.12) 


n-n c \ ( i* \ (4.22) 

.X = f y- j 

Similarly, the double sum is determined by equating Equa- 
tions (4.2) and (4.16): 

n~n c i (4.23) 

£ = 
i=l y=l 


H 

6 



10 


<i>„-2 = <!>(/« - 27) = 4, - 27’^ l > + 4 <£> - 


(4.13) 


T 3 

6 


7* 

24 


15 


$j*-3=<fr On ~3T) = 4>/t 


T 2 

2 


(4.14) 


27 ^~ <j) ( rt 3) + 81 74 


24 


20 


where T is the separation between intervals (“update inter- 
val”). By substituting Equations (4,10) through (4.14) in 
Equations (4.5) through (4.8), one can easily show the phase 
differences are related to input phase derivatives by 25 


A$„ tl =rC + -^- «?+-£- <t>< 3, +-§- €+■■ ■ 

(4.15) 


(4.16) 

> 

II 

r 

+ 

(4.17) 

A>„ 1 = r'C ) + ... 

(4.18) 


By respectively equating Equations (4.15) through (4.18) 35 
with Equations (4.1) through (4.4), one obtains of 
unknowns, where the unknowns are and the loop a set 
of equations whose number is equal to the number sums. 
Thus, these unknowns can be expressed in terms of the 
derivatives of input phase. 

To complete initialization of the loop, an estimate of 
starting model phase must be computed. To be exact, this 
estimate must account for steady-state phase error. For an 
arctangent phase extractor, tracking error is equal to residual 
phase (neglecting system-noise error and possible cycle 
ambiguities). For a sine phase extractor, however, steady- 45 
state tracking error becomes 


(<j)„-<}) m X j=arc sin (2 tc5 (J)^)/27C (4. 1 9) 

where phase is measured in cycles. Thus, model phase for 50 
the n th interval, after accounting for tracking error, is given 
by 


{ <j )„ - 5<j)jj for an arctan extractor (4.20) 

<J)„ - arcsin(27t5<|) J j)/27C for a sine extractor 55 

where phase is measured in cycles. 

The solution for a fourth-order loop will be presented in 
detail. If we assume that time derivatives of (j)(t) are negli- 
gible above <j)„ (4) , and equate Equations (4.18) and (4.4), we 
obtain a value for steady-state residual phase for a fourth- 60 
order loop: 



The loop sums are determined in a similar fashion. By 65 
equating Equations (4.3) and (4.17), one can solve for the 
single sum: 


k a 




n—n c 

Z 

i=l 



where the single sum is known from Equation (4.22). 
Finally, the triple sum is determined by equating Equations 
(4.1) and (4.15): 


(4.24) 


TL 

24 


n-n c n~n c j 

I £ £ 8% 

1=1 r=l /=! 


where the single and double sums are determined by means 
of Equations (4.22) and (4.23). Thus, if a priori values of 
signal phase and its derivatives are available at the start of 
a track, all of the loop variables can be initialized through 
use of Equations (4.20) through (4.24). If the a priori 
information is sufficiently accurate, the loop will start track- 
ing inlock, with no transients. 

Based on these equations for loop sums and model phase, 
the loop can be initialized as follows. Suppose a priori 
estimates of signal phase and its derivatives for a phantom 
“n* A interval” at time t„ are available. Then Equations (4.22) 
through (4.24) can provide estimates for loop sums at 
interval “completion” while Equation (4.20) provides an 
estimate of model phase at interval center, with steady- state 
phase error given by Equation (4.21). To project ahead to the 
first interval to be processed, the (n-tT) r/ * first estimate phase 
rate according to Equation (4.1) using these estimates for 
loop sums and steady-state phase error. Model phase for the 
(n+l) //z interval is then computed according to Equation 
(2.2) on the basis of this estimated (n+l)** phase rate and n !h 
a priori signal phase. The resulting values for the (n+l) /h 
phase and rate are then used to initialize the NCO for the 
(n+l)' A interval in the usual fashion. At the completion of the 
(n+1)'* interval, sums are updated in the standard fashion 
and substituted in Equation (3.1) of Equation (3.2) to 
provide the (n+2)^ estimate of rate. (If n c =l, the estimated 
steady-state phase error 5(j>^ is used as residual phase a 
second time). Equation (2.2) then provides an estimate of the 
(n+2) /; * model phase. Loop iteration is normal after this 
point, and steady-state lock should be achieved. 

The above analysis has focused on a fourth-order loop. 
Similar analysis can be carried out for loops of order lower 
than fourth, and the results are presented in Table 4-1. The 
preceding derivation assumed a signal with polynomial 
phase so that a steady-state phase error would develop. 
Under less ideal dynamics, the above initialization process 
will not eliminate transients but can greatly assist direct 
acquisition with higher-order loops. Similarly, if the deriva- 
tives of signal phase are known, but phase is not, the loops 
sums can be initialized using Equations (4.22) to (4.24), with 
initial loop phase arbitrarily set to zero. Again, loop acqui- 
sition will be greatly enhanced. 



5,602,883 


17 


18 


5.0 TWO MEASURES OF LOOP 
PERFORMANCE 

5.1 MEAN TIME TO FIRST CYCLE SLIP 

5 

Simulations have been carried out to determine mean time 
to first cycle slip, (T Uf ), for loops with phase and phase-rate 
feedback, no computation delay, supercritical damping, and 
a sine phase extractor with perfect amplitude normalization. 

The tracking-error criterion for detecting a cycle slip was 10 
l(j)-<j)J>0.75 cycles. After each slip, the loop was reinitial- 
ized with perfect initial values so that it would start off in 
steady-state lock with no transients. A Gaussian random- 
number generator simulated noise for the counter-rotation 
sums. 15 

Loops of 2 nd to 4 th order have been simulated on the basis 
of the loop constants in Table 3-5. Assumed values of B/T 
ranged between 0.02 and 2.0 and loop SNR between 0 and 
6 dB. Results are summarized in FIGS. 12a and 12 b, where 
B^ times mean time to first cycle slip is plotted versus loop 20 
SNR. In terms of cycle slips, loop performance deteriorates 
somewhat as loop order increases, given a fixed loop SNR. 

For a given loop order and loop SNR, however, cycle-slip 
performance improves as B L T increases, as shown in FIGS. 

12a and 12 b where B^(T U ,) is plotted versus B £ T for 2 nd to 25 
4 th order loops, given a loop SNR of 7 dB. With a 3 rd order 
loop, for example, FIGS. 13a and 13b indicates that B^T^,) 
improves by two orders of magnitude when B L T is increased 
from 0.02 to 0.5. 

As a test of the simulation software, the cycle-slip criteria 30 
have been made the same as Viterbi’s criteria in his exact 
closed form solution (5) for a first-order loop in the CU limit. 

To within a statistical error of about 10%, our results for 
BUTi^) agree with Viterbi’s theoretical predictions up to 
SNR^M- dB, the maximum loop SNR attempted. 35 

0 5.2 STEADY-STATE PHASE ERROR 
Loop performance at large values of B L T has also been 
assessed in terms of the steady-state phase error (SSPE). In 
the CU limit, SSPE is proportional to B^ for an N //l order 
loop, as can be derived from Tables 4-1 and 3-2. For large 40 
values of B^T, however, the SSPE in a root-controlled loop 
does not decrease as B L ~ N as B L increases, as illustrated in 
FIGS. 13a, 13b, 14 a and 146. FIGS. 13a, 136, 14a and 14 b 
plot as a function of B L T the dimensionless coefficient 
required to multiply the CU-limit form for SSPE. These 45 
plots pertain to loops of order 1 to 4, with phase and 
phase-rate feedback, with supercritical damping or standard 
underdamping, and with the indicated computation delay. At 
B L T=0, the coefficient is equal to the CU-limit value. As B/T 
increases, the increase in th is coefficient relative to the 50 
zero-B/T value is a measure of the “excess” SSPE relative 
to the nominal CU-limit values. As FIGS. 13a and 13 b 
indicate, for example, the SSPE at B/IM3.5 for a third-order, 
standard underdamped loop, is about four times as large as 
the CU limit would predict. 55 

FIGS. 13a, 13 b, 14 a and 14 b plot effective loop band- 
width as determined from SSPE, where “effective” denotes 
the decrease in bandwidth relative to the B^ model. In 
FIGS. 13a and 13 b, for example, the effective bandwidth is 
about 0.6 times the actual loop bandwidth when B/T^.5 for 60 
a second-order loop. 

SUMMARY 

A first-principles analysis of DPLLs has led to a new 65 
approach for parameterizing loops. Loop constants are com- 
puted on the basis of loop noise bandwidth and new con- 


trolled-root independent parameters that have direct physi- 
cal significance relative to decay times and damping. In the 
continuous-update limit, loop constants of loops of first to 
fourth order are obtained in closed form as a function of 
these new parameters. In a solution for a discrete-update 
(DU) loop, however, complexity of the equations leads to a 
numerical approach. The analysis has been applied to loops 
with phase and phase-rate feedback, given either zero com- 
putation delay or a computation delay equal to one update 
interval. With the new parameterization, exact selection of 
loop bandwidth and damping behavior can be carried out for 
high order loops, even when B/T is large. 

Simulations of loop behavior in terms of mean time to first 
cycle slip have been carried out for loops of first to fourth 
order based on the new parameterization, including the 
larger values of B^T. For a given loop bandwidth, loops with 
larger B L T exhibited a considerably better (larger) mean 
time to first cycle slip than those with smaller B L T values. 
Loop performance has been assessed on the basis of steady- 
state phase error (SSPE). As B/T increases or a given value 
of B^, SSPE is essentially constant or small B L T values (e.g., 
B^T^O.02) but increases for larger values of B L T. Plots 
provide a measure of this “excess” SSPE for large B L T 
values for loops with phase and phase-rate feedback. 


APPENDIX A. LOOP BANDWIDTH FROM 
LOOP CONSTANTS IN THE 
CONTINUOUS-UPDATE LIMIT 

We wish to solve for loop bandwidth 

f ^ 

2 B l = dv\H{v)\ 2 

in terms of loop constants K f , where the transfer function is 
given by 


Kh^TW)*- 1 + K2(i27W) yv " 2 + . . . + K N (A.2) 

H(V) ~ (i2-kv) n + Ki^TCV)^ 1 + k 2 (j7jw)^ 2 + ... + k n 

Under the change of variable u=2 tcv, we have 


Ki + K2(iuy^ 2 + . . . + Kjv 
{iuf + + Kniiu'f- 2 + . . . + k n 

(A.3) 

r do r °° 

4tc5z, = J du\H(u)\ 2 = J duH*(u)H(u) 

(A.4) 

From Gradshteyn & Ryzhik (3.112) (I. S. Gradshteyn and 
I. M. Ryzhik, Table of Integrals, Series and Products, 
Corrected and Enlarged Edition, Academic Press, New York, 
P. 218 (1980), we have 

f 00 ^ 8N(u) i7Z M n 

J U «0 A; V 

(A.5) 

where 


g N (u) - bjU 2(N ^> + b 2 u 2(N ~ 2 > + ... + b N 
h N (u) = a d^ + a jU N ~ J + . . . + a N 

(A.6) 

(A.7) 

and 
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with a computation delay of one update interval and phase 
( A -8) and phase-rate feedback becomes 

n - 1 n - 1 i (B.l) 

A<j>m,n+1 + K\§ mj x-\ + Ki Z <jVi + K-i Z Z <j) fflt / + • • • = 

5 z=l i= 1 7=1 

n-1 n-1 * 

+ K 2 Z (j)j + Ks Z Z <j)j + . , . 
i=l i=l 7=1 


Based on this equation, we find in a fashion similar to 
10 Sections 3.3.1 through 3.3.3 that the closed-loop transfer 
function is given by 


r 7M D(z)-z(z-D N W 

H(Z) ~ ~D(z) 

15 where the denominator is defined by 

Z)(z^zk-i) N +(z--i)^ 1 ^ 1 +za~i) w - 2 J fi: 2 +z 2 a-i) w - 3 A:3+ ...+z N ~ 
i K n (B.3) 


Comparing Equations (A.3) and (A.4) with (A.5), we get an 2Q 
expression for g^u): 

gu(u) = [Ki(iu) N ~ l + K 2(iuf~ 2 + . . . + Ktf] ( A - 10 ) 

[Ki(-m)^* 1 + k i(-iuf- 2 + . . . + Ktf] 

N r min(N~U~l) 1 

= Z tflHl K, 2 + 2 r {-l?Kj- k Kj +k 
7=1 L &=1 J 

and an expression for h^(u): 

M“)=(+^w^(+0^ 1 K/" 1 +(+i) w ' 2 K2t/ v " 2 + . . .+k n (A. 11) 30 


where the choice of -l-i rather than -i has been made to place 
the roots of h^u) in the upper half-plane. Finally, by 
equating like terms in (A.6) with (A. 10) 
and (A.7) with (A. 11), we have the coefficients a j and b,- in 35 
terms of the loop constants icy. 


(+if, if j - 0; 

otherwise. 


bj = TCj 2 + 2 


Z (-\) k Kj^Kj +k 

h= 1 


(A..12) 


(A. 13) 


40 


one can now express loop bandwidth as a function of loop 
filter constants by combining (A.4) and (A.5) to obtain 

45 

i Mn (A. 14) 

4ao A/y 

where the right side is given by Equations (A. 8) and (A.9) 
with the substitution of Equations (A. 12), and (A. 13). 

Table 3-1 lists results of Equation (A. 14) and presents B L 50 
as a function of for the first four loop orders. The table 
also presents B L as a function of the a’s after substitution of 
Equation (3.12). 


55 

APPENDIX B. ADAPTATION OF ANALYSIS TO 
LOOPS WITH COMPUTATION DELAY OF ONE 
UPDATE INTERVAL 


This polynomial has N+l roots in contrast to the zero- 
computation-delay case in Equation (3.42) and can be rep- 
resented by 


D(zMz-zi)(z-z 2 )(z-z 3 ) ■ ■ • (z-Zam-iH) (B.4) 

Upon expansion of D(z), a constraint equation for the extra 
(N+l)" root comes from the coefficient of z N when Equa- 
tions (B.3) and (B.4) are equated: 

N (B.5) 

ZN+i=N-Zzi 
1 

Thus if roots z x , z 2 , . . . , z N are designated as selectable 
roots, the (N+l)" dependent root is determined by this 
equation. As explained in Section 3.3.6, we are free to place 
the N selectable roots, as before, according to loop noise 
bandwidth, damping parameter and relative-magnitude 
parameter. The DU analysis of a computation-delay loop is 
otherwise parallel to the analysis for zero-computation- 
delay loop. The linear equations relating the loop constants 
to the independent roots are similar to Equations (3.49) 
through (3.53) and will not be shown here. Results for this 
implementation have been presented in the text. 

TABLE 3-1 


Loop Bandwidth from Loop Constants in the 
Continuous-update Limit 


l sl order 


Kl 

Bz = — 
2 nd order 


B/. = - 


Ki 2 + K 2 

^T“ 


3 rd order 


(1 + a 2 ) 


Bl — 


Ki 2 1 C 2 - K1K3 + k 2 2 

4 (kjK 2 - k 3 ) 


)- 4 - 


a 2 - a 3 + a 2 2 
a 2 -ct3 


4 th order 


The preceding DU analysis is based on a loop filter with 60 
no computation delay. Even though the analogous derivation 
for a loop filter with a computation delay of one update 
interval is in many ways a straightforward generalization, 
there are a few differences that are worth mentioning. In 
analogy with Equation (3.34), the loop equation for a loop 


Ki 2 K 2 K3 - K1K3 2 - Ki 3 K4 + K 2 2 K3 - K;K 2 K4 “ K3K4 
4(KlK 2 K3 - K 3 2 - Ki 2 K4) 

ki a 2 ot 3 - OC 3 2 - 04 + a 2 2 ct 3 - a 2 ot 4 - <*304 
4 a 2 a 3 - (X3 2 - 0C4 
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Loop Filter Constants in the Continuous-update Limit 


1 st order 
Kj = 4B l T 
2 nd order 

Ki = 4B l T- 


1 


1 + (X 2 
Controlled-root; 
Standard: 

3 rd order 

a 2 -(X3 


K 2 = a 2 K!2 

1 — Tli 2 


«2 = 


1 


K] = 4B^T - 


a 2 - a 3 + ax 

Controlled-root: 

Standard: 

4 th order 


K 2 = a 2 Ki 2 

2X 2 + a-m 2 ) 


0t2 = - 


K 3 = a 3 K, 3 

x. 2 a -ni 2 ) 


(2 + Xi) 2 


a 3 = - 


(2 + X 2 ) 3 


-4 


Ki = 4B l T- 


a 2 a3 - a 3 z - 04 


a 2 «3 - a 3 2 - 04 - a 2 C4 ~ 


K 2 = a 2 Ki 2 


K3 = a 3 Ki 3 


(*304 + ara 3 


Controlled-root: 

Standard: 


« 2 = - 


4X 2 + (1-71i 2 ) + 

^2 2 (1 ~ T12 2 ) 

(2 + 2 X . 2 ) 2 


«2- — 


a 3 = - 


«3 = 


2X 2 (1-Th 2 ) + 

2?t 2 2 (l-T1 2 2 ) 

(2+2X2? 


K4 = ouKi 4 


04 = 


A. 2 2 (l - tu 2 )(1 -TI2 2 ) 


(2 + 2 X 2 ) 4 




TABLE 3-3 


Loop-filter Constrants for Typical Implementations 
in the Continuous-update Approximation 


Loop constants Conventional parameters 


K 2 K 3 K 4 r k a 

K, (OjK, 2 ) (a 3 K, 3 ) (a 4 K, 4 ) (l/otj (<x+) (a„r 3 ) 


Supercritically damped: ti, 2 - 0, \ = 1, for all roots 


1 st Order 

4 B l T 





2 nd Order 

-5 

t k > 2 

4 



3 rd Order 

32 

y Ki 2 

4-k, 3 3 

1 

3 


4 th Order 


T k ' 2 

1 If 3 1 Tf 4 8 

16 Kl 256 Kl 3 

4 

9 

2 

27 


Standard underdamped: T\? = -1, ^ = 1, for all roots 



1 st Order 

4 B l T 





2 nd Order 

T BiT 

T K > 2 

2 



3 rd Order 

-§-B t T 

jK, 2 

2 rf 3 9 

27 Kl 4 

3 

8 
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TABLE 3-3-continued 



Loop-filter Constants for Typical Implementations 
in the Continuous -update Approximation 




Loop constants 


Conventional parameters 

K, 

K 2 

(c^K, 2 ) 

k 3 

(a 3 K! 3 ) 

k 4 

(a 4 K t 4 ) 

(1/Oj) 

k 

(a 3 r 2 ) 

a 

(a 4 r 3 ) 

4 th Order 64 

27 

B t T -j- K, 2 

t k ' 3 

-^ K > 4 

2 

1 

2 

1 

8 


TABLE 3-4 

Loop Bandwidth from Loop Constants for Discrete- update Loops with 
Phase and Phase-Rate Feedback and No Compoutation Delay 

1 st order 

K] 

BlT ~ 2(2 -Ki) 

2 nd order 

2K 1 2 + 2K 2 + K 1 K 2 
BlT ~ 2Ki(4-2Kj-K 2 ) 

3 rd order 

4Ki 2 K 2 - 4KiK 3 + 4K 2 2 + 2K!K 2 2 + 4K] 2 K 3 + 4K 2 K 3 + 3KiK 2 K 3 + K 3 2 + K]K 3 2 
Bz - T " 2 (KiK 2 - K 3 + KiK 3 )(8 - 4Ki ~2k 2 - K 3 ) 

4 th order 


8 Ki 2 K 2 K 3 - 8 KiK 3 2 - 8 KAK 4 + 8 K 2 2 K 3 - 8KiK 2 K4 - 8K 3 K4 + 4K!K 2 2 K 3 + 8 K! 2 K 3 2 + 8 K 2 K 3 2 + 4K! 2 K 2 K4 + 8 K 2 2 K 4 - 
20KiK 3 K4 - 8 K 4 2 + 6 K]K 2 K 3 2 + 2K 3 3 + 14Kt 2 K 3 K4 + 14K 2 K 3 K4 + 4 KiK 2 2 K4 - IOK 1 K 4 2 + 1 1KiK 2 K 3 K4 + 
2KiK 3 3 + 5 K 3 2 K 4 + 7 K! 2 K 4 2 + 6 K 2 K 4 2 + SK^IQ + 5 KiK 2 K 4 2 + 4 K 3 K 4 2 + 4 KiK 3 K 4 2 + K 4 3 + K 1 K 4 3 
2 (KjK 2 K 3 - K 3 2 - Ki 2 K4 + KiK 3 2 + KiK 2 K4 - 2K 3 K4 + 2K1 1K 3 K4 - 
K 4 2 + KiK4 2 )(16 - 8 K 1 - 4K 2 - 2K 3 - K 4 ) 


TABLE 3-5 


Loop Filter Constants for a Discrete-update Loop with Phase 
and Phase-rate Feedback and Supercritically- damped Response 


b l t 

1st order 
Ki 

2nd order 


3rd order 


4th order 


Ki 

k 2 

Ki 

k 2 

k 3 

Ki 

k 2 

k 3 

k 4 





No computation delay 





0.001 

0.003992 

0.003193 

2.553e-06 

0.002903 

2.812e-06 

9.084e-10 

0.002747 

2.833e-06 

1.299e-09 

2.234e-13 

0.005 

0.01980 

0.01582 

6.309e-05 

0.01438 

6.941e-05 

1.11 8e-07 

0.01361 

6.989e-05 

1.597e-07 

1.369e-10 

0.01 

0.03922 

0.03130 

0.0002488 

1 0.02845 

0.0002733 

1 8.778e-07 

0.02692 

0.000275 

1.251e-06 

2.138e-09 

0.02 

0.07692 

0.06125 

0.0009677 

0.05567 

0.001060 

6.765e-06 

0.05269 

0.001065 

9.617e-06 

3.265e-08 

0.03 

0.1132 

0.08993 

0.002118 

0.08174 

0.002312 

2.220e-05 

0.07738 

0.002321 

0.0000312 

1.578e-07 

0.05 

0.1818 

0.1438 

0.005576 

0.1307 

0.00605 

9.485e-05 

0.1237 

0.006059 

0.0001337 

1.113e-06 

0.075 

0.2609 

0.2051 

0.01176 

0.1864 

0.01267 

0.0002936 

0.1766 

0.01265 

0.0004113 

5.055e-06 

0.1 

0.3333 

0.2607 

0.01965 

0.2369 

0.02101 

0.0006405 

0.2245 

0.02094 

0.0008915 

l.439e-05 

0.15 

0.4615 

0.3572 

0.03931 

0.3248 

0.04147 

0.001847 

0.3080 

0.04113 

0.002540 

5.978e-05 

0.2 

0.5714 

0.4379 

0.06264 

0.3983 

0.06523 

0.00378 

0.3780 

0.06443 

0.005139 

0.0001570 

0.25 

0.6667 

0.5061 

0.08835 

0.4606 

0.09089 

0.006432 

0.4375 

0.08942 

0.008652 

0.0003222 

0.3 

0.7500 

0.5643 

0.1155 

0.5139 

0.1175 

0.00976 

0.4885 

0.1152 

0.01300 

0.0005671 

0.35 

0.8235 

0.6142 

0.1436 

0.5598 

0.1444 

0.01371 

0.5327 

0.1411 

0.01808 

0.0008996 

0.4 

0.8889 

0.6575 

0.1720 

0.5998 

0.1712 

0.01821 

0.5712 

0.1668 

0.02380 

0.001325 

0.45 

0.9474 

0.6952 

0.2006 

0.6348 

0.1977 

0.02321 

0.6050 

0.1920 

0.03006 

0.001844 

0.5 

1.0 

0.7282 

0.2291 

0.6657 

0.2235 

0.02864 

0.6349 

0.2166 

0.03679 

0.002459 

0.6 


0.7828 

0.2851 

0.7173 

0.2732 

0.04059 

0.6852 

0.2634 

0.05133 

0.003967 

0.7 


0.8257 

0.3394 

0.7585 

0.3196 

0.05371 

0.7258 

0.3069 

0.06692 

0.005832 

0.8 


0.8599 

0.3915 

0.7920 

0.3629 

0.06769 

0.7590 

0.3471 

0.08318 

0.008026 

0.9 


0.8874 

0.4414 

0.8196 

0.4030 

0.0823 

0.7865 

0.3842 

0.09983 

0.01052 

1. 


0.9096 

0.4890 

0.8426 

0.4402 

0.09735 

0.8097 

0.4184 

0.1167 

0.01328 

1.2 


0,9425 

0.5779 

0.8782 

0.5065 

0.1283 

0.8462 

0.4790 

0.1503 

0.01951 

1.4 


0.9645 

0.6588 

0.9042 

0.5636 

0.1597 

0.8734 

0.5308 

0.1832 

0.02650 
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TABLE 3-5-conlinued 





Loop Filter Constants for a Discrete-update Loop with Phase 
and Phase-rate Feedback and Supercritically- damped Response 




1st order 

2nd order 


3rd order 



4th order 


b l t 

K, 

K x 

k 2 


k 2 

k 3 

Ki 

k 2 

k 3 

k 4 

1.6 


0.9793 

0.7328 

0.9237 

0.6130 

0.191 

0.8942 

0.5754 

0.2151 

0.03410 

1.8 


0.9889 

0.8006 

0.9386 

0.6559 

0.222 

0.9106 

0.6141 

0.2457 

0.04216 

2. 


0.9950 

0.8631 

0.9502 

0.6935 

0.2525 

0.9236 

0.6479 

0.2749 

0.05058 

2.5 


1.0 

1.0 

0.9697 

0.7689 

0.3258 

0.9466 

0.7158 

0.3419 

0.07270 

3. 




0.9811 

0.8248 

0.3947 

0.9612 

0.7667 

0.4011 

0.09566 

3.5 




0.9881 

0.8672 

0.4591 

0.9710 

0.8058 

0.4533 

0.1189 

4. 




0.9925 

0.8997 

0.5194 

0.9778 

0.8366 

0.4996 

0.1422 

4.5 




0.9953 

0.9248 

0.576 

0.9827 

0.8613 

0.5409 

0.1652 

5. 




0.9971 

0.9444 

0.6291 

0.9864 

0.8814 

0.5779 

0.1879 





Computation delay = 1 

update interval 




0.001 

0.003976 

0.003181 

2.538e-08 

0.002892 

2.796e-06 

9.015e-10 

0.002737 

2.816e-06 

1.288e-09 

2.211e-13 

0.005 

0.01942 

0.01554 

6.135e-05 

0.01414 

6.749e-05 

1.077e-07 

0.01339 

6.795e-05 

l,537e-07 

1.306e-10 

0.01 

0.03779 

0.03023 

0.0002357 

0.02752 

0.0002588 

8.165e-07 

0.02606 

0.0002604 

1.163e-06 

1.952e-09 

0.02 

0.07177 

0.05734 

0.0008737 

0.05224 

0.0009552 

5.897e-06 

0.04951 

0.0009599 

8.364e-06 

2.745e-08 

0.03 

0.1027 

0.08185 

0.001832 

0.07461 

0.001993 

1.809e-05 

0.07076 

0.002000 

2.554e-05 

1.231e-07 

0.05 

0.1571 

0.1245 

0.004476 

0.1136 

0.00482 

7.032e-05 

0.1078 

0.004819 

9.836e-05 

7.611e-07 

0.075 

0.2148 

0.1685 

0.008742 

0.1538 

0.009282 

0.0001955 

0.1462 

0.009237 

0.0002700 

3.007e-06 

0.1 


0.2046 

0.01371 

0.1869 

0.01433 

0.0003895 

0.1778 

0.01420 

0.0005309 

7.609e-06 

0.15 


0.2594 

0.02487 

0.2377 

0.02515 

0.0009724 

0.2265 

0.02467 

0.001289 

2.610e-05 

0.2 




0.2740 

0.03595 

0.001784 

0.2618 

0.03495 

0.002295 

5.906e-05 

0.25 




0.3000 

0.04617 

0.002793 

0.2878 

0.04452 

0.003473 

0.0001076 

0.3 







0.3071 

0.05316 

0.004765 

0.0001722 

0.35 







0.3211 

0.06074 

0.006129 

0.0002543 


TABLE 3-6 




Loop Filter Constants for a Discrete-update Loop with Phase 






and Phase-rate Feedback and Standard-underdamped Response 




1st order 

2nd order 


3rd order 



4th order 


b l t 

K, 

K, 

k 2 

Ki 

k 2 

k 3 

Kx 

k 2 

k 3 

k 4 





No computation delay 





0.001 

0.003992 

0.002661 

3.545e-06 

0.002603 

3.016e-06 

1.311e-09 

0.002367 

2.804e-06 

1.661e-09 

4.923e-13 

0.005 

0.01980 

0.01319 

8.752e-05 

0.01290 

7.437e-05 

1.611e-07 

0.01173 

6.907e-05 

2.037e-07 

3.008e-10 

0.01 

0.03922 

0.02609 

0.0003448 

0.02552 

0.0002926 

1.264e-06 

0.02321 

0.0002716 

1.594e-06 

4.693e-09 

0.02 

0.07692 

0.05106 

0.001338 

0.04996 

0.001133 

9.720e-06 

0.04545 

0.001051 

1.222e-05 

7.147e-08 

0.03 

0.1132 

0.07499 

0.002922 

0.07338 

0.002469 

3.156e-05 

0.06679 

0.002289 

3.955e-05 

3.447e-07 

0.05 

0.1818 

0.1199 

0.007658 

0.1174 

0.000445 

0.0001355 

0.1070 

0.005966 

0.0001687 

2.420e-06 

0.075 

0.2609 

0.1713 

0.01607 

0.1677 

0.01345 

0.0004177 

0.1530 

0.01243 

0.0005159 

1.093e-05 

0.1 

0.3333 

0.2179 

0.02670 

0.2133 

0.02226 

0.0009073 

0.1949 

0.02054 

0.001112 

3.096e-05 

0.15 

0.4815 

0.2991 

0.05288 

0.2929 

0.04370 

0.002595 

0.2683 

0.04022 

0.003136 

0.0001273 

0.2 

0.5714 

0.3675 

0.08345 

0.3599 

0.06842 

0.005269 

0.3305 

0.06282 

0.006285 

0.0003313 

0.25 

0.6667 

0.4258 

0.1166 

0.4171 

0.09491 

0.008901 

0.3838 

0.08697 

0.01049 

0.0006738 

0.3 

0.7500 

0.4760 

0.1511 

0.4662 

0.1222 

0.01341 

0.4299 

0.1118 

0.01562 

0.001176 

0.35 

0.8235 

0.5196 

0.1862 

9.5089 

0.1496 

0.01871 

0.4701 

0.1367 

0.02156 

0.001851 

0.4 

0.8889 

0.5577 

0.2213 

0.5463 

0.1768 

0.02470 

0.5056 

0.1613 

0.02817 

0.002705 

0.45 

0.9474 

0.5914 

0.2561 

0.5793 

0.2034 

0.03129 

0.5370 

0.1853 

0.03534 

0.003738 

0.5 


0.6214 

0.2902 

0.6085 

0.2294 

0.03838 

0.5650 

0.2087 

0.04296 

0.00495 

0.6 


0.6721 

0.3560 

0.6590 

0.2788 

0.05381 

0.6128 

0.2532 

0.05920 

0.007884 

0.7 


0.7134 

0.4181 

0.6983 

0.3247 

0.07047 

0.6521 

0.2944 

0.07632 

0.01145 

0.8 


0.7475 

0.4763 

0.7316 

0.3671 

0.08796 

0.6848 

0.3325 

0.0939 

0.01558 

0.9 


0.7762 

0.5306 

0.7593 

0.4061 

0.1060 

0.7124 

0.3675 

0.1116 

0.02021 

1. 


0.8007 

0.5813 

0.7829 

0.4421 

0.1243 

0.7360 

0.3997 

0.1293 

0.02527 

1.2 


0.8399 

0.6727 

0.8205 

0.5059 

0.1612 

0.7743 

0.4569 

0.1639 

0.03644 

1.4 


0.8699 

0.7524 

0.8490 

0.5603 

0.1979 

0.8038 

0.5058 

0.1971 

0.04869 

1.6 


0.8936 

0.8223 

0.8713 

0.6072 

0.2337 

0.8272 

0.5479 

0.2286 

0.06171 

1.8 


0.9127 

0.8839 

0.8891 

0.6478 

0.2684 

0.8461 

0.5846 

0.2582 

0.07525 

2. 


0.9286 

0.9385 

0.9035 

0.6832 

0.3019 

0.8618 

0.6167 

0.2860 

0.08913 

2.5 


0.9587 

1.051 

0.9297 

0.7544 

0.3802 

0.8910 

0.6817 

0.3482 

0.1245 

3. 


0.9827 

1.134 

0.9470 

0.8076 

0.4510 

0.9111 

0.7310 

0.4013 

0.1599 

3.5 




0.9592 

0.8485 

0.5151 

0.9256 

0.7695 

0.4470 

0.1945 

4. 




0.9680 

0.8805 

0.5735 

0.9366 

0.8004 

0.4967 

0.2281 

4.5 




0.9747 

0.9061 

0.6269 

0.9452 

0.8256 

0.5213 

0.2605 

5. 




0.9798 

0.9267 

0.6760 

0.9520 

0.8465 

0.5519 

0.2917 
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TABLE 3-6-continued 





Loop Filter Constants for a Discrete-update Loop with Phase 
and Phase-rate Feedback and Standard-underdamped Response 




1st order 

2nd order 


3rd order 



4th order 


b l t 

Kj 

K, 

k 2 

K, 

k 2 

k 3 

K, 

k 2 

k 3 

k 4 





Computation delay = 1 update interval 




0.001 

0.003976 

0.00265 

3.516e-06 

0.002594 

2,997e-06 

1.299e-09 

0.002358 

2.784e-06 

1.645e-09 

4.859e-13 

0.005 

0.01942 

0.01299 

8.487e-05 

0.01270 

7.229e-05 

1.549e-07 

0.01155 

6.717e-05 

1.959e-07 

2.862e-10 

0.01 

0.03779 

0.02533 

0.0003248 

0.02474 

0.0002769 

1.170e-06 

0.02252 

0.0002573 

1.478e-06 

4.259e-09 

0.02 

0.07177 

0.04827 

0.001194 

0.04709 

0.001019 

8.391e-06 

0.04293 

0.0009475 

1.058e-05 

5.935e-08 

0.03 

0.1027 

0.06919 

0.002479 

0.06730 

0.002121 

2.554e-05 

0.06164 

0.001971 

3.210d-05 

2.636e-07 

0.05 

0.1571 

0.1060 

0.005937 

0.1030 

0.005094 

9.769e-05 

0.09425 

0.004732 

0.0001220 

1.594e-06 

0.075 

0.2148 

0.1448 

0.01130 

0.1401 

0.009726 

0.0002657 

0.1285 

0.009020 

0.0003292 

6.120e-06 

0.1 


0.1775 

0.01725 

0.1709 

0.01489 

0.000518 

0.1571 

0.01379 

0.0006363 

1.503e-05 

0.15 


0.2300 

0.02963 

0.2193 

0.02572 

0.001236 

0.2021 

0.02372 

0.001494 

4.860d-05 

0.2 


0.2713 

0.04155 

0.2554 

0.03630 

0.002167 

0.2360 

0.03335 

0.002579 

0.0001037 

0.25 


0.3071 

0.05222 

0.2833 

0.04616 

0.003232 

0.2622 

0.04227 

0.003795 

0.0001781 

0.3 




0.3055 

0.05518 

0.004373 

0.2832 

0.05037 

0.005075 

0.0002685 

0.35 




0.3242 

0.06338 

0.005530 

0.3002 

0.05766 

0.006374 

0.0003710 

0.4 







0.3143 

0.06420 

0.007664 

0.0004821 

0.45 







0.3263 

0.07006 

0.008927 

0.0005986 
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TABLE 4-1 


Transient-free Initialization of DPLLs Based on APriori Information 
For a loop filter specified by update interval T, computation 
delay n c update intervals, constants K lt K 2 , . . . , K N and input 
signal phase known in the form <j>(t) = + (t - tj + ... + 

— ■ <j)^(t - t n ) N at time t„, loop sums and loop model phase (j) m?n 
are initialized at the end of the “phantom” n* interval as indicated. 


1 st Order 


= 

2 ad Order 


T 

Ki 


U) 


til 


.2 5<fc 
i=l 

3 rd Order 


T~ 

K 2 






n-n c 

i=l 


Z 2 8^ 

i=l j=l 


4 5h Order 


T 3 

K 3 




(T 2 ^ - Ki&Jijj) 




n-n c 

Ki&])„-Ka £ 
i=l 


&§ss 


X 5 kj>,- 

i=l 


-p4 

"kT 




~T 


49- 


Ka&fc 


n-n c { 
i=l j=l 


- j ^- 1 T 2 ^” 1 — IT ” 


n-n c 


1-lie ! i 

£ Z 1 
i=l j=l k=l 


For all orders: 

f i>»- 


1 

K4 




T 2 


>-^3 rp 4 n— n c n— n c j 

-4- 't> < n 3, + -4r- C’-K^-K, Z &t>,-K 3 S Z&fo 

0 1=1 1=1 J=1 


4> rt - Arcsin(27t b$ ss )/2 k. 


Atan extractor 
Sin extractor 


Having thus disclosed a preferred mode of carrying out 50 
the invention, what is claimed is: 

1. A method for implementing a digital phase-locked loop 
having discrete updates to loop parameters and a pre- 
selected loop order greater than 1 and characterized by 55 
preselected parameters for loop bandwidth, computation 
delay, update-interval damping and decay rate; the method 
comprising the steps of: 

a) determining the roots for the characteristic equation of 60 
the loop from the preselected parameters; 

b) determining the loop constants from the loop roots of 
step a); and 

c) implementing a digital loop filter incorporating the loop 65 
constants of step b) thereby meeting the requirements 
imposed by said preselected loop parameters. 


2. The method recited in claim 1 wherein said digital 
phase-locked loop is of the type having immediate update 
for modifying feedback for the very next update interval. 

3. The method recited in claim 1 wherein said digital 
phase-locked loop is of the type having delayed updates for 
modifying feedback for a subsequent update interval after 
the very next update interval. 

4. A method for establishing the loop filter constants in a 
digital phase-locked loop of order greater than 1, with 
discrete updates to loop parameters and characterized by 
preselected parameters for loop bandwidth, computation 
delay, update-interval, damping and decay rate, the loop of 
the type having a numerically-controlled oscillator for gen- 
erating counter-rotation phasors based upon phase-related 
feedback for complex multiplication of a digitized input 
signal, the loop of the type also having an accumulator of 
counter-rotated signal phasors and a residual phase extrac- 
tor; the loop filter of the type estimating next phase rate and 
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generating feedback for input to the numerically-controlled 
oscillator; the method comprising the steps of: 

a) determining the roots for the characteristic equation of 
the loop from the preselected parameters; 

b) determining the loop constants from the loop roots of 
step a); and 

c) implementing a digital loop filter incorporating the loop 
constants of said loop filter as determined in step b) 
thereby meeting the requirements imposed by said 
preselected loop parameters. 

5. The method recited in claim 4 wherein said digital 
phase-locked loop is of the type having immediate update 
for modifying feedback for the very next update interval. 

6. The method recited in claim 4 wherein said digital 
phase-locked loop is of the type having delayed updates for 
modifying feedback for a subsequent update interval after 
the very next update interval. 

7. A method for implementing a digital phase-locked loop 
having discrete updates to loop parameters and a pre- 
selected loop order of 1 and characterized by preselected 
parameters for loop bandwidth, computation delay and 
update-interval; the method comprising the steps of: 

a) determining the root for the characteristic equation of 
the loop from the preselected parameters; 

b) determining the loop constant from the loop root of step 
a); and 
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c) implementing a digital loop filter incorporating the loop 
constant of step b) thereby meeting the requirements 
imposed by said preselected loop parameters. 

8. A method for establishing the loop filter constant in a 
5 digital phase-locked loop of order 1 with discrete updates to 
loop parameter and characterized by preselected parameters 
for loop bandwidth, computation delay and update-interval, 
the loop of the type having a numerically-controlled oscil- 
lator for generating counter-rotation phasors based upon 
10 phase-related feedback for complex multiplication of a 
digitized input signal, the loop of the type also having an 
accumulator of counter-rotated signal phasors and a residual 
phase extractor; the loop filter of the type estimating next 
phase rate and generating feedback for input to the numeri- 
15 cally-controlled oscillator; the method comprising the steps 
of: 

a) determining the root for the characteristic equation of 
the loop from the preselected parameters; 

20 b) determining the loop constant from the loop root of step 
a); and 

c) implementing a digital loop filter incorporating the loop 
constant of said loop filter as determined in step b) 
thereby meeting the requirements imposed by said 
25 preselected loop parameters. 

* # & * * 



